Natural Inflation, Planck Scale Physics 
and Oscillating Primordial Spectrum 



Xiulian Wang"^, Bo Feng", Mingzhe Li", Xue-Lei Chen''''^, Xinmin Zhang" 
" Institute of High Energy Physics, Chinese Academy of Science, P.O. Box 918-4, Beijing 100039, P.R. China 
National Astronomical Observatories, Chinese Academy of Sciences, 

20A Datun Road, Beijing 100012, P. R. China and 
" Kavh Institute for Theoretical Physics UCSB, CA 93106, USA 

In the "natural inflation" model, the inflaton potential is periodic. We show that Planck scale 
physics may induce corrections to the inflaton potential, which is also periodic with a greater fre- 
quency. Such high frequency corrections produce oscillating features in the primordial fluctuation 
power spectrum, which are not entirely excluded by the current observations and may be detectable 
- - - in high precision data of cosmic microwave background (CMB) anisotropy and large scale structure 

I (LSS) observations. 
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I. INTRODUCTION 



In the past decade, inflation theory has successfully passed several non-trivial tests. In particular, recent cosmic 
microwave background (CMB) observations show that the spatial geometry of the observable Universe is very close 
, to flat just as inflation theory predicts. Inflation theory also offers an elegant way of generating the 

^ ' primordial fluctuations which seed the formation of galaxies and large scale structures (LSS) (see e.g., Ref. Q and 
I references therein). In particular, slow-roll inflation models predict that the perturbations are adiabatic, gaussian, 
^ ' and nearly scale-invariant (i.e with a power index Us ~ 1). It has long been known that these predictions are also in 
fSJ ' broad agreement with the observed properties of large scale structures and CMB anisotropy, although at present the 
. data is still not very restrictive ,7]: = O.OIIq Jj. 

This general success of inflation theory brings up the hope of extracting even more detailed information of the 
, inflaton potentially,^] from high quality observational data. There are two complementary approaches to this problem. 
' In the first approach, one tries a model-independent (subject to some standard assumptions) reconstruction of the 
] primordial power spectrum, then the inflaton potential from the data [oL llCt ITll| . Alternatively, one can look for 

■ specific features in the power spectrum and study their observational consequences. In particular, there have been 

■ many investigations on inflationary models with broken scale invariance ll^ 1 1 7t ITsj l . Such features have 
been invoked to explain the tentatively observed feature at A: ~ 0.05AfDc~^p^ l2(il2lLl23| . or even to solve the small 



Oh; 

O ■ scale problem of the CDM model |2^.^ 



In the present work, we consider a new type of feature, periodic in the primordial power spectrum. This type of 



^ ' feature is interesting from both a theoretical and phenomenological point of view. Theoretically, if discovered, it gives 
. . [ a strong hint on the nature of the inflaton field. Phenomenologically, it might change the position, shape or even the 
^ number of acoustic peaks in the CMB power spectrum. 

k>( _ This paper is organized as follows: in Sec. II, we describe our model, and show how this type of feature could arise 
. from Planck scale physics by constructing a toy model. Our toy model, which is based on the "natural inflation" 
' model, is by no means the only possibility, but in this context it is particularly easy to see how this might happen. 
In Sec. Ill we derive the power spectrum in this model, and then consider how it would affect CMB and large scale 
structure in Sec. IV. The final section, Scc.V is on the summary and discussions of our results. 



II. THE MODEL 



In addition to phenomenological success, a compelling infiation model should also be based on plausible particle 
physics. To be in agreement with observations, the infiaton potential must be very flat. Since the radiative corrections 
to a scalar field mass are quadratically divergent, some physical mechanism, e.g., a symmetry is required to maintain 
the flatness of the potential, unless we want to accept ad hoc fine tuning. 

For the natural inflation model (S^ , several possible physical mechanisms 'sTI are available for producing a inflaton 



^ For other solutions to this problem, see Refs. |24| - |2S| and for a recent review on this issue see Ref. 1291 
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potential of the form 

y(0)=A4[l±cos((/.//)]. (1) 

Customarily, the positive sign is taken, with identical results for the negative sign. If A ~ 10^^ GeV and / ^ 10^^ GeV, 
which are the grand unification scale and Planck scale respectively, a successful inflation model could be obtained, 
with the correct quantum fluctuation amplitudes |22]- 

In this paper we extend the natural inflation model above and include a term in the potential in Eq.JQ) which is 
also periodic but with greater frequency: 

Vpianck (0) = 5K^ cos(iV0// + 9) . (2) 

We will argue below that this term may come from the Planck scale physics. 

In natural inflation, one introduces a pseudo Nambu-Goldstone boson (PNGB) as the inflaton. Specifically let's 
consider an anomalous global U{\) symmetry which is spontaneously broken at scale / via a non- vanishing value of a 
complex scalar field '-^ :^/^*"^^'^ with the Nambu-Goldstone boson. At lower energy scale A like Axion a potential 

for the Goldstone boson (p is generated from the non-perturbative effects and its mass is order of ~ 

Recently there are a lot of interests in studying the effects of the Planck scale physics on the primordial spectrum 
and CMB anisotropy |33j . In the original paper of this subject, Martin and Brandenberger considered a modified 
dispersion relation and showed explicitly a sizable effect. In this paper we take an effective lagrange approach to new 
physics and argue that Planck scale physics can correct the inflaton potential, generate oscillating and scale-dependent 
power spectrum which may be detectable in the future. 

The Planck scale physics which can be parameterized by higher dimensional operators makes two types of contribu- 
tions to the effective lagrangian of the inflaton cf). Since is the phase of the scalar field if the Planck scale physics 
preserves the U{\) symmetry the Goldstone theorem requires that there must be derivatives involved in the operators 

with (j). The operators with the lowest dimension are like . The authors of Ref. [s^ have studied the effects of 

pi 

this type of operators on CMB and shown that the corrections to CMB is order of Hf^ji^^-^^/Mpi^^^^/. ^ 10"-^^. And 

Ref. pointed out that higher order terms like ^^"^^4 '^^ would be possible to enhance the effects. In this paper we 

pi 

study another type of the Planck scale physics effects where the higher dimensional operators involve no derivatives 
and modify the inflaton potential. This requires an explicit violation of the global symmetry U(l). In fact, as argued 
in Refs.lSG] the quantum gravitational interactions break the global symmetries and one expects the existence of a 
type of higher dimension operators like 



(3) 

Mpi 



which gives rise to a correction to the inflaton potential of 0: 

Vpianck (0) = ^ A* COs{Ncj)/f + 0) , (4) 

where 

S = 2-'^-^/'+'\gMN\^U^) (5) 



A4 \Mpi^ 

with 9 the phase of the parameter gMN- Taking M — 2, we have 

where we have taken j^l ^ 1. 

For most of the studies in this paper we will assume that the gravitational interaction correction to the mass of the 
Goldstone scalar is less than that from the non-perturbative effects i.e. SN^ < 1. This implies that N in general is 
expected to be large. For instance taking / = O.bMpi and A = 2.6 x 10~^Mp;, SN^ < 1 requires that N > 36, and for 
/ = 0.95Mp;, A = 9.7 X 10~'^Mpi , one needs N > 90. This requirement makes the basic feature of the original natural 
inflation phenomenologically unchanged, however theoretically it requires an convincing argument to understand why 
the operators with dimensions less than N are forbidden or much suppressed. In this paper we will focus on the 
phenomenological studies. 
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The inflaton potential under investigation in this paper is the sum of Vpianck and the one in Eq. |^ , namely 

= A4[1 + cos(0//) + <5cos(7V0// + 0)]. (7) 

In additional to the two parameters A and / in the original natural inflation we introduce in our model three more 
parameters 6,N and 6. The S and N, however are related by Eqs.(5) and (6). And for simplicity of the discussions 
we will restrict ourself to two cases that = or 6* = tt in this paper. Choosing 6* = tt is equivalent to changing S 
to ~S. Before concluding this section, we point out that S term shifts the minimum of the potential away from zero, 
however since S is extremely small, it will not change the results when minimizing the potential to zero, which we 
have made a numerical check for the calculation of the primordial spectrum and CMB results. 



III. PRIMORDIAL POWER SPECTRUM CALCULATION 

During inflation, the contribution of other matter can be neglected, and the background evolution of the Universe 
is described by 



(8) 



^ + 3ff0 + ^0 = 0, (9) 
where H = a/a is the expansion rate, and = dV/ dip. Slow rolling (SR) requires 

(10) 

~ H^^ IGtt^V ' Stt V ^ ^ 

For the model under consideration in this paper, we have 

(sin(0//) + 57Vsin(A^0//))2 



(12) 



167r/2 (l + cos(0//))' 

^^16./H l + cos(<^//) )■ 



The primordial spectrum of the scalar perturbation is defined as 

Psik)^^Kk\\ (14) 



where C,k is the coefficient of the Fourier transform of the Bardeen parameter C, |37l |. The scalar spectrum index is 
given by 

dlnPJk) , , 

= l + (15) 



The primordial spectrum of the tensor perturbation is 



fc3 



PT{k)^j:^\M\ (16) 

where ipk is the coefficient of the Fourier transform of tensor linear perturbations ip US ■ ^'^"^ tensor spectrum 
index: 

dln/Mfc) 
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In SR regime, 

ns = l-4e-2/3, tit = -2e. (18) 

One can see from Eq. H12|l and Eq. that the relative contributions of the 6 term to the SR parameter e given by 
the terms oi 6N or are proportional to 1/N ot 5 for a fixed 5N'^. Numerically they are small with the parameters 

(large N and small S) we consider in this paper. But (3 can be strongly modulated by large 6N'^, then Ps{k) and Ug 
are changed correspondingly. Note here V and are changed minorly by S term, hence the background evolution of 
the inflaton field will not be affected much. 

During inflation, the equation of motion for a Fourier mode of fluctuation is [s^ 

< + (^fc" - ^ j uk - 0, (19) 

where 

u = zC, z= a^/H. (20) 

Here the prime denotes the derivative with respect to the conformal time d-q = dt/a{t). Well inside the horizon, 
according to flat spacetime quantum field theory, the vacuum modes are 

Wfc-^^e"'*^" as fc2^— ; (21) 



/2k z 
Coming out of the horizon, the growing mode solution is 

z" 

Uk (X z for fc^ < — , (22) 

z 

with no explicit dependence on the behavior of scale factor a in this limit (frozen). For the slow roll approximation 
(SRA), ^ can be expressed in terms of two SR parameters e and f3: 

^ = 2a»H»(l + |9 + , + i^^ + l,^+ll, + llffl. (23) 
If e and (3 can be approximatively considered as constants, Eq. (|19|l can be rewritten as 

< + - ^~^) = ^^^^ 
where = + \ and i] ^ fn (i^) ' "^^^ solution of Eq. is a Bessel function, 

Uk = iv^e'('^+^)^(-r,)^i/(i)(-fcr;). (25) 
When k/aH 0, from the asymptotic form of H^, one obtains an expression of the scalar spectrum Ps: 

where Ps{k, a) varies with time. Traditionally one takes k = aH, and gets the primordial spectrum 



4^2 r2(|) 



Psik) -^^^Al- ^r-'^ \k=aH . (27) 



The calculation for the gravitational wave spectrum is very similar |87l |. The equation of motion for tensor linear 
perturbations is 

< + {k' - ^)vk = 0, (28) 
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where v = aip, and 



Vk'^^e''''^ as (29) 



Vk (xa for < — . (30) 

a 



Therefore, 



For SRA, one can get 



where 



^T(fc) = ^|^P. (31) 



/'T(fc)-2^-^^(l-.)^-^^U..., (32) 



^ (33) 



1 -e 2 

The Stewart-Lyth formula given in Eq. (|27|l and Eq. H32() are vahd under certain conditions that both of the parame- 
ters e and P are almost constants and the asymptotic value of Ps{k) can be approximated by Ps{k = aH). In Ref . j39|| 
the authors pointed out a subtlety in the approximations in obtaining the formula in Ea. (|27|l . however numerically 
they have also demonstrated that the Stewart-Lyth formula works fairly well for the chaotic and natural inflation 
models. In the presence of the S term, the SRA is not necessarily true and (3 is oscillating, especially in the case when 
SN'^ is relatively large. Thus, to ensure the validity of Ea. H18|) . we calculate the spectrum numerically in this paper. 

We numerically solve Eq. ((HJ, Eq. 0, Eq. 119|l and Eq. (|28(l and obtain the horizon crossing amplitude for each 
mode with the following initial conditions 39, 40] : 

..(0) ^ -1= , (34) 



4(0) = ^ ; (35) 



..(0) ^ -= , (36) 



-1.(0) = ^ . (37) 



2k 

Our results show that for SN'^ < 1, the Stewart-Lyth formula works well and the relative error is within 1% in 
comparison with the numerical results. However, for 6N'^ > 1, which we will discuss in detail in section V, the error 
could be large. In Fig^ we plot Psik) with numerical results and the Stewart-Lyth approximation. One can see that 
the difference is about 30 percent; for gravitational waves, the difference is little. 

In the natural inflation model, the scalar spectrum index Ug is generally smaller than 1. The S term in our inflaton 
potential induces a modulation on the power spectrum. In Fig. 12 we plot the primordial scalar spectrum and index 
for a typical set of model parameters: / = 0.7Mpi, and A"' — 2.0 x 10~"'^^Mp;. As can be seen from the figures, the 
amplitude of the modulation in this model could be 0(10%). 

Before concluding this section we should point out that the size of the modulation depends on the parameters of 
the model, and Ps(fc) behaves quite differently for small and large /. This can be understood qualitatively from the 
following two equations for natural inflation in the SR regime (Here we've set SnG = 1) : 

isin(0//) ^^^^ 



din a 1 + cos((/)//) ' 
and 



^/2f _ ^r 

(f>j^ ^2f&YCsin(—^=J=e ^), (39) 
^1 + 2/2 



6 



where Af is the number of e- folding between the corresponding horizon exit and the end of inflation . For smaller 
/, 4>_\f would be much smaller than for larger /. For example (j)fif « 2.03 for / — 0.95Mp;, while / = OAMpi gives 
(j>^ « 6.3 X 10"'' for a given number of e-f aiding M = 70. Such a small (j) for / — QAMpi would also make ^j^^ 
much smaller than / — 0.95Afp;. This is also true for our model when SRA is well satisfied. 

In Figl^l we plot the Ps(fc) and ns{k) for a small /. One can see from this figure that for / — QAMpi, Ps{k) has 
lost the feature of oscillation in the range of scale relevant to CMB and LSS, but the amplitudes of Ps and ns change 
a lot. 



IV. RESULTS AND COMPARISON WITH OBSERVATIONS 



In this section we study the implications of our model on CMB and LSS. The matter power spectrum is related to 
the primordial scalar power spectrum by 

P{k)^T\k)Ps{k), (40) 

where T{k) is the matter transfer function, which can be calculated for a given set of the cosmological background 
parameters. The matter power spectrum in a large range of scales can now be obtained by combining different data 
sets of CMB and LSS measurements 41]. For example, currently the Lyman alpha forest probes comoving scale of 
0.1 ~ lOhMpc-^, the 2dF galaxy correlation 0.01 OAhMpc^^, and CMB measurements 0.001 OAhMpc^^. The 
CMB anisotropy is related to the primordial scalar power spectrum bv.42l] 

(A(ni)A(n2)) = ^^(^^ + ^)CiPi[cos{ni ■ ^2)], (41) 

(=0 



where Ti{k) is the photon transfer function, which can be calculated for a given set of the cosmological parameters. 
Here Cj, C; and T; stand for the temperature angular power spectrum and transfer function. We've omitted the 
subscript T for simplicity. We calculate Ci by modifying the publicly available CMBFAST code. Similar 

calculations for running ns have been performed in, for example, Refs. 16, 32, 45j. The fiducial model adopted in 
our calculation is the best fit model of Ref. 0: /i = 0.64, = 0.66, nbh^ = 0.02 and Qk = 0. We used the combined 
CMB data set from Ref.^^ in our plots. In Fig. 01 we plot the CMB angular power spectrum and the matter power 
spectrum with the primordial power spectrum shown in Fig[21 One can see that the theoretical predictions of our 
model are quite different from that with S — Q and the CMB data disfavor the case with SN'^ = 0.8. 

One interesting point of our model is the running behavior of rig in the case with a small /, which we have shown 
in Fig. El One can see that without S term the model predicts below 0.8, which is ruled out already by the data. 
However with the help of the S term the theoretical predictions can be made to be consistent with the observations 
which can be seen from Fig.[5l Numerically in natural inflation with / — OAMpi it predicts ns ~ 0.76. And the fine 
tuning of the Vq would not be able to make the theory consistent with C; and P{k) 7]. However the presence of the 
d term makes the fit better. 

Up to now we have restricted ourself to the parameter space SN'^ < 1, however it would be interesting to take the 
potential in Eq. ((7|) as a phenomenological model and study the cosmological implications with an extended parameter 
space SN"^ > 1 . In this case, there are four parameters, where we take / , N and 6 free, however A can be normalized 
by observations. 

The effect of the modulation on the power spectrum depends on the amplitude, frequency, and phase. And the 
effect is more apparent for a high frequency (large N), as can be seen in Figs. \6lWi In Fig.jSlwe take N = 599 and 
find the variation of primordial scalar spectrum index is very frequent and its amplitude changes from about 0.3 to 
1.6, although the value of S is very small. However, our results also show that the tensor perturbation Pt changes 
little since e <C 1. We have checked the validity of the consistency relationj48j^. One can see from Fig. [3 for different 
range of fc, {nxl is larger or smaller than Pt/Ps- The reason responsible for the violation of consistency relation 



The consistency relation holds for single-field SRA inflation. The authors of Ref. |4S| argue that the trans-Planckian physics may change 
the vacuum and provides an example for the violation of consistency relation. 
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in our model is due to a large running of /3 from —0.3 to 0.3 , which indicates SRA isn't satisfied well. However, 
we will show below that this type of primordial spectrum doesn't contradict to the current observations. In Fig|S| 
we plot the CMB anisotropy and the matter power spectrum for the given parameters and find that there can be a 
clear modulation on the matter power spectrum, which is somewhat similar to the baryonic oscillation wiggles, but 
has an entirely different origin. Of course, at present there is no observational evidence of such wiggles in the power 
spectrum. The peaks found in recent CMB measurements seem to agree reasonably well with the predictions of power 
law primordial spectrum. In the large scale structure data, there is some tentative report on wiggles in the power 
spectrum ,50^,.5i,_52, 53]. These are generally thought to be due to baryonic oscillation, but the effect discussed here 
could also give rise to such wiggles. Since the position and frequency of the baryonic wiggles can be predicted, and 
there is no reason to expect that the primordial power spectrum modulation to coincide with it, there is hope to 
distinguish these two cases with high precision data. However, we note that at present time there is still no clear 
evidence even for the baryonic wiggles [s^. l55|. Note that the modulation on the power spectrum is coherent, with 
different Ps the shape of the CMB peaks can be also quite different, some would drastically change the structure of 
the Ci spectrum, the peaks of the Ci can be split into two in some particular cases. In the left panel of FiglHl the 
first peaks have been clearly split, while in FigEI it shows that the original first and second peaks have been split, 
the third peak is highly raised and new peaks are also generated. In Fig. ^| we plot the CMB polarizations which are 
given by 



"""'WTTiI t^'V-'n'V-'ir'sik) (14) 

for the E-mode and cross correlation spectrum polarization. These effects might be detectable by further CMB 
observations [5^. 

Current CMB observations do not exclude entirely the oscillating primordial spectrums. Consider the difference 
between C° and (7/ where (7° is calculated with / = 0.95Mpi, 5 = (ns « 0.945) and Cl with / = 0.95Mpi , N = 599 
, S = —3 X 10^^, then compare them with the cosmic variance [STj limits on C". In Fig^] the solid line stands for 

C^-C° ■ ... .... ~i ~n 

' ' and the region between the dashed lines is given by the cosmic variance limits. Here C; and Cj are normalized 
by COBE. Cosmic variance plays a fundamental limit to the C; spectrum that can be measured, with 



One can see from Fig^Jthat Cf and Cf are hardly distinguishable. Even though the primordial spectrum (see, FiglS)) 
is oscillating, its contribution to CMB, because of being averaged over multiple 1, becomes insignificant in this case. 



V. DISCUSSIONS AND SUMMARY 



In this paper we present a model which is a variation of natural inflation. We have shown two features of the 
primordial spectrum of this model, oscillating and scale-dependence and studied the implications on CMB and LSS. 
In the presence of the S term the parameter space allowed for a successful natural inflation will be enlarged relative 
to the original natural inflation modeljs^. When the parameter space is enlarged and extended to SN^ > 1 for 
example due to some other physical motivations _5gJ, there are several additional interesting effects. Although the 
SR approximation is violated and the spectral index oscillates with a large scale variation, there could be a large 
parameter space not ruled out by the observations. As we can see from Fig|Slwhen we gradually increase the value of 
6 the effects on CMB firstly take place on the first peak which can be slightly split, meanwhile wiggles on the matter 
power spectrum are gradually enhanced. The effects on CMB first peak and large scale structure are potentially 
observable and can be tested by future precise experiments. 

In summary we have studied in this paper a model which is a variation of natural inflation and show that there are 
some interesting phenomenological features of this model, such as oscillating and scale-dependence in the primordial 
spectrums. And we have also discussed their implications on CMB and LSS. 
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FIG. 3: The same as Fig|21 however with different model parameters, / — OAMpi, A* — 5.0 x 10 ^"^ Mpi. From above, the 
hnes represent (5iV^=-0.7, -0.2, -0.07, 0, 0.07, 0.2, 0.7 respectively. 




FIG. 4: The CMB anisotropy and matter power spectrum for the parameters shown in Fig. |21 with 5N'^ = 
0.8, 0.42, 0, -0.42, -0.8 from above. The observational LSS data is from the PSCZ catalogue^. 
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FIG. 5: Ci (left) and matter power spectrum (right) for a fixed / = OAMpi but different values of otfier model parameters. 
For the solid line, A'' = 1.3 x IQ-^^M^,, SN'^ = 0; The dotted line stands for A'' = 3.6 x lO-'^^Mj.i, SN^ = -0.34; For the 
dashed line, A"* = 5.2 x W'^^Mpi, SN'^ = -0.78. 




In[k/kj 

FIG. 6: The primordial power spectrum index for / = 0.95Mpi, N = 599, 5 = -3 x 10 
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FIG. 7: Plot of Inrl and Pt/Ps for / = 0.95Mp,, N = 599 and 5 = 



-3 X 10"^ 
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FIG. 8: CMB anisotropy and matter power spectrum for / = 0.95Mp;, A'' = 599. Left panel: the dash dot dotted, the dashed, 
the dotted, the thicker sohd and the thmncr sohd Une stand for 5 = —3 x 10^'', 5 — 3 x 10^', 5 — —3 x 10^^, 5 = and 
(5 = — 8 X 10~^ respectively; Right panel: The dash dotted dot, the thicker solid and the thinner solid stand for 5 = —3x 10""', 
(5 = -8 X 10"^ and 5 = -3 x 10"^ respectively. 




FIG. 9: Theoretical prediction on Ci with / = .95Mpi, 5 = 0.0007 and N = 799. 




FIG. 10; CMB polarization for the same parameters as in the right panel of Fig. 
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FIG. 11: Plot of AC; and the cosmic variance limits. 



